Project supported by the National Natural Science Foundation of China (Grant Nos. 51677008, 51377182, 51707028, and 11647098), the Fundamental Research Funds of the Central Universities, China (Grant No. 106112017CDJQJ158834), and the State Key Development Program for Basic Research of China (Grant No. 2014CB541602).
Project supported by the National Natural Science Foundation of China (Grant Nos. 51677008, 51377182, 51707028, and 11647098), the Fundamental Research Funds of the Central Universities, China (Grant No. 106112017CDJQJ158834), and the State Key Development Program for Basic Research of China (Grant No. 2014CB541602).
† Corresponding author. E-mail:
Project supported by the National Natural Science Foundation of China (Grant Nos. 51677008, 51377182, 51707028, and 11647098), the Fundamental Research Funds of the Central Universities, China (Grant No. 106112017CDJQJ158834), and the State Key Development Program for Basic Research of China (Grant No. 2014CB541602).
The conventional magnetic resonance imaging (MRI) equipment cannot measure large volume samples nondestructively in the engineering site for its heavy weight and closed structure. In order to realize the mobile MRI, this study focuses on the design of gradient coil of unilateral magnet. The unilateral MRI system is used to image the local area above the magnet. The current density distribution of the gradient coil cannot be used as a series of superconducting nuclear magnetic resonance gradient coils, because the region of interest (ROI) and the wiring area of the unilateral magnet are both cylindrical side arc surfaces. Therefore, the equivalent magnetic dipole method is used to design the gradient coil, and the algorithm is improved for the special case of the wiring area and the ROI, so the X and Y gradient coils are designed. Finally, a flexible printed circuit board (PCB) is used to fabricate the gradient coil, and the magnetic field distribution of the ROI is measured by a Gauss meter, and the measured results match with the simulation results. The gradient linearities of x and y coils are 2.82% and 3.56%, respectively, less than 5% of the commercial gradient coil requirement.
Conventional magnetic resonance imaging (MRI) must work in highly homogeneous field,[1,2] but the volume and weight of the equipment restrict its applications, especially for huge object. For conventional closed magnet structures, the sample must be placed inside the magnet. However, the unilateral magnet provides a new direction of application for MRI.
Unilateral magnetic resonance imaging (UMRI) is a new type of nuclear magnetic resonance (NMR) measurement method. We use the characteristics of unilateral magnets to build a system to realize magnetic resonance imaging (MRI), which is called UMRI. Compared with conventional closed MRI, the UMRI can be easily moved for its open structure, small size, and the sample can be nondestructively tested from the surface above sample, obviously UMRI has a good prospect.[3,4] In recent years, the research on imaging systems based on UMRI equipment has also made some progress.[5–7] Blümich presented a two-dimensional (2D) phase-encoding imaging method to realize 2D imaging with a unilateral NMR probe.
The gradient coils play an irreplaceable role in MRI, which should be specially designed for unilateral magnet. At the same time, since the unilateral magnet has a static gradient field perpendicular to the surface of the magnet (the z-axis direction as shown in Fig.
The design of the gradient coil is an inverse problem of the electromagnetic field. The design method can be divided into regular separation winding method and distributed winding method.[8,9] The distributed winding method does not predetermine the shape of the coil, but rather calculates the desired current density distribution by a predetermined magnetic field distribution within the region of interest (ROI), and then use conductive copper plates or distributed winding to simulate the current density distribution, in order to determine the specific shape and size of the coil winding. This method is also easy to determine inductors, energy consumption and self-shielding constraints and other conditions. The design method includes stream function method,[10,11] target field method,[12,13] harmonic coefficient method,[14] and equivalent magnetic dipole method.[15,16] Owing to different magnet structures, the gradient coil design method is slightly different. At present, some design methods are applied to the design of coils in Halbach magnets,[17,18] which has some similarities to the unilateral magnet gradient coil design. However, in this paper the design of gradient coils for unilateral magnets is still lacking. We adopt the equivalent magnetic dipole method to design gradient coils, focusing on the characteristics of a unilateral magnet, and whose ROI and wiring area are both cylindrical side arc surfaces.
In this study, a unilateral MRI magnet was designed for mobile UMRI. The main magnetic field is parallel to the upper surface of the magnet. The default direction is the X direction as shown in Fig.
In this paper, the wiring area of the gradient coil is a cylindrical side arc surface, and it can be divided into q units, each unit is a tiny cylindrical piece with thickness h, length a, and arc length a, as shown in Fig.
According to the law of current continuity, important preconditions for current densities can be obtained as
Therefore, we can define a hypothetical scalar stream function S(
The stream function defines the function with the continuous direction changing in a limited area. Therefore, the stream function is suitable for the current distribution problem discussed in this paper. The stream function is used to describe the current density distribution of the base current in each split cylindrical surface, wherein the contour distribution of the stream function can represent the actual current routing, and the difference in stream function between the adjacent contour lines is the actual current value
We can solve the value of the stream function S(
Therefore, substituting Eq. (
In Fig.
If a ≪ |
According to Eq. (
Therefore, the relationship between the target magnetic field
If the design of the coil only makes the magnetic field in the ROI meet the target magnetic field, the structure of the gradient coil may be not smooth, which will increase the resistance of the coil and make the local temperature too high. Hence, the energy loss needs to be optimized. The energy consumption expression produced by all magnetic dipole elements can be written as
The optimal mathematical model is constructed as follows:
The constraints in Eq. (
The latter item on the right-hand side in Eq. (
After obtaining the optimal value of the stream function Sq, the coil structure can be obtained by drawing the contour of the stream function. In addition, the optimal value of the penalty parameter λ needs to be determined. This requires the gradient linearity δ to be used as an optimization indicator, which is expressed as
(i) Set the initial value of λ to be λ = 1010.
(ii) Obtain the optimal solution Sq from Eq. (
(iii) Calculate the Bx value by substituting Sq into Eq. (
(iv) Set the optimal exit conditions for δ < ε, where ε is the gradient linearity requirement. When exit conditions are met, the loop ends. The current value λ is used as the optimal penalty parameter, or λk + 1 = aλk, where a is usually set to be 10[i]. Then repeat steps (ii)–(iv).
According to the actual magnet structure and simulation parameter optimization, the wiring area of the gradient coil in this paper is a 200-mm-length cylindrical side arc surface. Its radius is 45 mm, and its central angle is 120°. Wiring area is divided into 50 subareas equally in the X, Y direction, and is divided into a total of 2500 units, each unit can be approximated as a square, split unit side length a total. The ROI is a cylindrical surface with a radius of 35 mm central angle of 90°, length of 60 mm. We divide the cylindrical surface into 100 area elements, each corresponding to a target field point. Using the above analysis we can carry out the calculation to obtain the X and Y gradient coil.
When designing the gradient coil, at λ = 1010, we set the parameter ε = 3 to exit from the loop to obtain X, Y gradient coil, and the structure of the gradient coil can be obtained as shown in Figs.
As the conventional coil resistance is too large, more serious heat is produced, we use the anti-coil structure to reduce the coil resistance. The final structure is shown in Fig.
According to the simulated coil structure, the anti-coil structure design and optimization, and finally using the flexible printed circuit board (PCB) processing, the gradient coil is obtained as shown in Fig.
This study uses the United States FW.Bell8030 Gauss meter, 3D stepper platform, Agilent6653A current source constitutes a measurement platform shown in Fig.
Then the results measured from X gradient coil and Y gradient coil and simulation results are analyzed and compared in Tables
We inject direct currents (DCs) of 1 A, 2 A, and 3 A into the X and Y gradient coil, and obtain the distribution of the central axial magnetic field of ROI, which is shown in Fig.
First of all, in this work we design a unilateral magnet in the field of medical local imaging of mobile MRI equipment and realize the arc-shaped target magnetic field distribution with a natural gradient in the area of 10 mm cylindrical side arc surface above the magnet.
Then based on the principle of the unilateral specific direction of the main magnetic field, the gradient coil design adopts the equivalent magnetic dipole method. According to the requirements for this unilateral magnet, the gradient coil wiring area and the ROI are designed to be cylindrical side arc surface, in this paper, the design method for the above situation has been improved and optimized.
Finally, we apply the coil design results to real coils, use mature measurement platform to measure the gradient coil, and find that the magnetic field distribution in the ROI is close to the simulation results. However, because of wiring problems such as asymmetry, there are still some deviations.
The mobile UMRI will be widely used, especially for the open imaging of large volume samples (such as localized breast imaging). Gradient coil design for unilateral magnets is an important research foundation for the study of the mobile UMRI, so it needs to be further studied. In the future work, we will further investigate the sequence and imaging results of UMRI in combination with the existing UMRI equipment.
[1] | |
[2] | |
[3] | |
[4] | |
[5] | |
[6] | |
[7] | |
[8] | |
[9] | |
[10] | |
[11] | |
[12] | |
[13] | |
[14] | |
[15] | |
[16] | |
[17] | |
[18] |